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In quantum physics, measurement error and disturbance were first naively thought to be simply 
constrained by the Heisenberg uncertainty relation. Later, more rigorous analysis showed that the 
error and disturbance satisfy more subtle inequalities. Several versions of universally valid error- 
disturbance relations (EDR) have already been obtained and experimentally verihed in the regimes 
where naive applications of the Heisenberg uncertainty relation failed. However, these EDRs were 
formulated for discrete measurements. In this paper, we consider continuous measurement processes 
and obtain new EDR inequalities in the Fourier space: in terms of the power spectra of the system 
and probe variables. By applying our EDRs to a linear optomechanical system, we conhrm that a 
tradeoff relation between error and disturbance leads to the existence of an optimal strength of the 
disturbance in a joint measurement. Interestingly, even with this optimal case, the inequality of the 
new EDR is not saturated because of doublely existing standard quantum limits in the inequality. 

PACS numbers: 03.65.Ta, 42.50.-p, 42.50.Lc 


Introduction — The uncertainty principle is one of the 
most fundamental features of quantum physics; it pre¬ 
vents us from measuring non-commuting observables si¬ 
multaneously with arbitrary accuracy. The original for¬ 
mulation by Heisenberg [ij was based on a thought ex¬ 
periment that measures the position of an electron with a 
7 -ray microscope. Nowadays the Heisenberg uncertainty 
relation (HUR) is known as 

caOb ^ \Cab\ ■ (1) 

where Cab = /2, ca = c{A) is the error in the 

measurement of the observable A of a quantum system 
and r]B = oi.B) is the disturbance exerted back onto a 
system observable B, which does not commute with A. 

Long after the original proposal of the HUR, more 
rigorous analyses of the measurement process incorpo¬ 
rating the measuring device, or probe, revealed that the 
HUR, when used naively as an error-disturbance relation 
(EDR), can be violated, e.g. 0. Ozawa proved a revised, 
universally valid EDR [3|, 1^ : 

CAVB + cactb + ctaVb > \Cab\ , (2) 

where a'\ and are variances of the observables A 
and B, respectively. Soon after the proposal of the 
Ozawa’s EDR inequali^, it has been generalized to the 
error-tradeoff relation [a, 1^, in which the disturbance r]B 
is replaced with cb and the inequality is symmetrized 
about the variables A and B. The Ozawa’s EDR can 
be tightened by using stronger geometrical inequalities 
7 -Is. The strongest one has been derived by Branciard 
a l9| and is given by 

^"a^b + ^aVb + ‘^\J^\^^B ~ ^AB ^aVb > Cab ■ (3) 


violated while the universal EDRs hold. The above EDRs 
were derived for measurements performed at instants of 
time, i.e., discrete measurements, and cannot be applied 
directly to continuous measurements: in a discrete mea¬ 
surement, one can define the error and disturbance by 
comparing quantities before and after the measurement. 
However, notions of’’before and after” are irrelevant to a 
continuous measurement. In this paper, we redefine error 
and disturbance as the differences between ’’without and 
with the measurement” and extend the above EDRs to 
the continuous measurements. 

Fundamental variables — In general, we consider the 
quantum measurement of a system achieved via coupling 
to a device, or a probe; the system and probe have sepa¬ 
rate Hilbert spaces. The projection postulate of quantum 
mechanics is applied directly to the probe — and acts 
indirectly on the system. We suppose that a system ob¬ 
servable A is the target of measurement, with the meter 
observable of a probe M that is projectively measured, 
and consider the disturbance back onto another system 
observable B. 

For a continuous measurement process, we define the 
error and disturbance operators with time-dependent 
Heisenberg operators as N{t) = M{t) — Ao{t) and D{t) = 
B{t) — Bo{t), respectively, where the subscript ”0” stands 
for free evolution, for which the system and probe are de¬ 
coupled. These are motivated by Ozawa’s original defi¬ 
nitions for discrete measurement d, , but differ in their 
interpretations. Specifically, for a discrete measurement, 
N and D compare variables between before and after 
the measurement. In the continuous case, our N and 
D compare the difference of variables between with and 
without measurements. The error and disturbance in the 
measurement of Aq are 


These EDRs have been experimentally tested with 


the three-state method 


and with the weak- 


measurement method I2-ll4ll and found that the HUR is 


e\it)^{N\t)) = {{M{t)-Aoit)r), (4) 

pUt)^{D^t)) = {{m-Bo{m- (5) 
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Here (■ ■ ■) denotes ensemble average or time average for 
a stationary system. The variances of observables Aq and 
Bq are 

a^(t)^({ioW-(ioW)}"), (6) 

(7) 

These definitions of error, disturbance, and variances are 
general, but here we treat a stationary measurement, 
where any two-point correlation function depends only 
on time difference. Then the quantities in Eqs. 
become independent of time. 

Error-disturbance relation in Fourier space — We 
first note that Eq. m continues to hold if we re¬ 
place error and disturbance quantities, ca and rjg, by 
their time-dependent versions defined in Eqs. ®-© 
— and commutator Cab by the equal-time commuta¬ 
tor [Ao{t),Bo{t)]. This is because the algebraic relations 
that had lead to the original inequality continue to hold 
in this situation. However, at a steady state, the inequal¬ 
ity is the same for all t, and we only have one condition 
governing the total fluctuations. In order to obtain fur¬ 
ther constraints, we shall assume the noise processes to 
be stationary, and are therefore characterized by their 
spectra. 

For the later use, we introduce a filter function T{t), 
whose shape can be chosen arbitrarily; the filtered error 
and disturbance operators are 

Af{t)= f dt'T{t-t')N{t')=M{t)-Ao{t) , 

J —CO 

V{t)= f dt'T{t-t')D{t')=B{t)-Bo{t) , 

J —CO 

while Ao{t), Bo{t), M{t), and B{t) are defined in a similar 
fashion using the same filter. The new error e can be 
computed using 

pCO J(~) 

e\^{N\t)) = j^ —\T[nfs,{n). (8) 



FIG. 1. Comparison of the EDR lower bounds. The curves 
indicate the locations at which the Ozawa (red short-dashed), 
Branciard (blue solid) and Heisenberg (green) inequalities are 
saturated. Here we assume |x(H)| = 1. 

where Sri, SaA, and Sas are the power spectra of D, Aq, 
and Bq, respectively. 

Since the filtering does not change the Hilbert space of 
each observable and not affect the proof of Eq. @, the 
EDR in Eq. ([2]) still holds for the filtered quantities: 

+ <^AVI3 > \Cab\ ^ (12) 

with the commutator. Cab = *([-do(t),Ho(t)])/2. 

Substituting Eqs. (|H]), CT- dTTl) and choosing T to 
be narrowbanded within a very narrow frequency bin 
near H, the filtered version of the Ozawa’s inequality in 
Eq. (fT^ can be written as [l5| 

^s,{n)s^{n) + ^/Sr{n)s^r^{n) + ^ 5 <,^( h ), s ^( h ) 

>;r|(x(H))|. (13) 

where we have defined (using staionarity) 


Here H is angular frequency and we defined 

/ CO 

dT0(T)r(r)e*^^ , 

-CO 

with the Heaviside step function 0(r) and the (one-sided) 
power spectral density Se = Snn of N [l^. Note that for 
a stationary state the time dependence of ca disappears. 
Similar to Eq. ®, VBi O'An a 13 , are given by 

pCO IQ 

= l —imfSriin), (9) 

poo IQ 

<^A = l -^irmfs.Am, ( 10 ) 

poo IQ 

^B = J^ (11) 


[io(ll), Sl(H')] = 2TTih5{n - O')x(H) . (14) 

Since Eq. mjp plies to all frequencies, it generalizes 
Ozawa’s EDR [l5j into a spectral relation for our contin¬ 
uous measurement. 

Analogously, the generalization of the Branciard EDR 
to the continuous case reads [l5j |: 

+ s^A{^)Sn{n) 

+ 2^s4n)Sn{n) {5.,(h)^^,(h) - |(x(fi))|"} 

■ (15) 

Defin ing Se = S^/S^a, Sr, = Sn/S^B, and x = 
^(x) / \/ S^A Sob I '"’'0 obtain normalized Ozawa and Bran- 
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ciard inequalities: 


^Js,{n)Sr,{n) + + ^Jsr,{n) > |x(f^)| , (i6) 

5 ,( 0 ) + 5 ^( 0 ) + 2 y 5 ,( 0 ) 5 ^( 0 )(l-|x( 0 )|') > |x( 0 )|' 

(17) 

Let us briefly comment on the spectral expression of 
the Robertson’s inequality [1^, aAcrs > 5o])|, 

which is the basic uncertainty relation that holds for ar¬ 
bitrary system observables and has nothing to do with a 
measurement. Applying the filters in Eqs. (flUl) and (fTTl) 
(this is just a mathematical operation and not on actual 
data), we have the Robertson’s inequality in the Fourier 
space, 


To be more specific, let us consider the probe realized 
by a coherent light field with power Iq and frequency ujq ■ 
According to the derivation in [isj, the Qg and Fg are 
proportional to the phase and amplitude quadratures of 
the optical field. Then we have the canonical response 
function Rqf{^) = 1. Since we do not consider a feed¬ 
back depending on the position of a test mass, we have 
F = Fg. The power spectra are 


Sqo{^) 


hc^ 

SIgUJg ’ 


5F„(fl) = 


SIghiUg 


and the cross-correlation vanishes, SpoQoi^) = 0. 

Since Rgri^) = 1 for our optical probe, using 
Eqs. (I^ - (I^ . we have the error and disturbance in the 
Fourier space 


lx(^^)l < 1 ■ (18) 

If the system saturates the Robertson’s inequality, the 
third term in Eq. (El) vanishes. In other words, the Bran- 
ciard inequality becomes the tightest when |x(n)| = 1. 

On the other hand, comparing Eqs. du and © , we 
find that the HUR is given by keeping only the first term 
on the LHS of Eq. (fTBl) : 


N{n) = Qg{n) + R,,{n)Fg{n) , 

Din) = Rp,in)Fg{n). 

Thus the power spectra of the error and disturbance are 
5, = Sq„ + \R,,f Sfo + 2Re [Ro..SfoQo] , (24) 

Sr, = \Rp.\^SFo ■ (25) 


s,in)Spin) > Ix(0)|2 . (19) 

As we shall illustrate later, there are the case where the 
Ozawa’s inequality in Eq. (fT6l) is supported by the second 
and third terms on the LHS, thereby allowing the HUR 
to be violated. In Fig. [U three inequalities, Eqs. (US, 
ED, and dE]), are compared, setting |x(ll)l = 1- 
Application to a linear optomechanical system — We 
now apply our formalism to a linear optomechanical sys¬ 
tem, in which the position of a test mass is measured 
with an optical probe. Namely, we specify the system’s 
observables A and B to be the position and momentum 
operators of the system, x and p, and the meter variable 
M of the probe to be the phase quadrature of light Q, 
whose unit is the same as x. We assume that the in¬ 
teraction between the probe and the system is the von 
Neumann type, with Hamiltonian in the interaction pic¬ 
ture given by Hi{t) = —xg{t)Fgit), where the subscript 
”0” denote free evolution and Fg is a generalized force. 

For this linear system, the time evolution of any linear 
variable can be written in terms of its response to gener¬ 
alized forces. In Fourier space, if we define the response 
function 

2 /'°° 

RxYm = j^j dT0(T)C'xy(r)e*^" , (20) 

which satisfies i?^y(H) = Rxyi—n), then the dynamics 
of the system is described by [15|, 

x{n) = xg{n) + R,,in)Fin), ( 21 ) 

Pin) = pgin) + Rp,in)Fin), ( 22 ) 

Qin) = Qgin) + RQFin)xin). ( 23 ) 


We model the test mass as a harmonic oscillator cou¬ 
pled with environmental thermal bath. Free evolution 
of the position and momentum of the damped oscillator 
without thermal contributions are given by, e.g. 


m 


Xgit)=e—B2 
_ p-K™t/2 

-Xgit) 


^ , Pini . , 

Xini COS LUmi H-SmCJrnC 


rnuJm 

pgit) = [-mUJmXiniSrnUJmt + PiniCOSLUmt] 

mKr, 


where m, k^, and Wm are the mass, the decay rate, and 
the resonant frequency of the oscillator. The initial val¬ 
ues are fini = io(0) and Pi^i = Po(0)- The commuta¬ 
tor decays as but it is compensated by thermal 

fluctuations from the external heat bath. The two-time 
commutation relation between the position operators at 
different times, t and t' , are computed straightforwardly 
and then from Eq. (I20p . the response functions in the 
Fourier space are obtained. 


Rxxin) 


Rxpin) 


Rpxi^) 


_L__ 

H -j- iOui -f j‘2^ 

. W I e-^^ 

4cUyj.j ( H -j- -t- iKjnj^ 

-Rxpi ^), 


n ajrfi J 

(26) 

e*® 1 

n T j 


where W = -\/4wTTjr7t^_and 9 = arctan[—/Cm/(2wm)]. 
Then using the formula [l5| 


[Xo(H),yoi'(H')] = -2Trihdin-n') {Rxvin) - Ryxi^')} . 

(27) 
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FIG. 2. EDR applied to a damped harmonic oscillator. For the Branciard inequality, each curve represents Se (blue, dotted), 
Sn (red, dashed), and the LHS of Eq. (1171) (magenta, solid), respectively. Eor the Heisenberg inequality, the LHS of Eq. (1191) 
(green, dashed-dotted) is plotted. The horizontal line is the right-hand side (RHS) of the inequalities, |x(f2)|^ = 1. From the 
left to the right, a = louJoH'mc^uj'^) = CTopt/5, (Topt, and CTopt x 20, respectively, k is fixed to p = K/uJm = 0.3. 


and comparing with Eq. (HI, we have 

v(0) =_ _ 

^ {(11 + + ^C^/4}{(11 - u^m? + «^/4} ■ 

(28) 

The power spectrum of test-mass zero-point fluctua¬ 
tions are given by the fluctuation-dissipation theorem, 
Sxo = 2?ilm[i?3,a,(n)]. Using Eq. ([26l), 


m{(ll -b + K^/4}{(11 - + Km/4} ’ 

(29) 

= . (30) 


Since (xini) = 0 and (pini) = 0, we have SaA = Sxo and 

Sas ~ ^po ■ _ 

Eirst we see that the Robertson’s inequality in Eq. m 
is always saturated for a damped harmonic oscillator as 
explicitly verified with Eqs. (EHi-dSQl). Consequently, for 
such a minimum uncertainty state, the third term in the 
LHS of Eq. (ITSl) vanishes, and the inequality, generally 
keeping nonzero correlation term SpoQoi reads 

Q I c _ ^Qo 4~ ‘^\^xx\‘^Spg + 2 Re[Ra;3;>5'fQQq] 

' 2ftIm[J?,,] - • 

(31) 

If SpgQg remains real-valued, using the Braginsky’s in¬ 
equality Sor,SFn — I'S'foQoP ^ and after some algebra 
(derived in [l5|), we obtain 


Se + Sr, > \/2 + {2-V2){a^ + 62) -b (3 - 2V2)a%^ , 

(32) 

where a = Sp^Qg/h and b = Re[i?a;a;]/Im[i? 2 ; 3 ,]. While it 
is possible for Sf„q„ to be complex, this usually leads to 
the modification of the test mass’s dynamics, which we 
shall not consider. Erom this equation, we see that the 
best one can do is to have a = b = 0, that is, S'fqQo = ^ 
and Re[i? 3 , 2 ,] = 0, which takes place on resonance, where 
the inequality is a factor of from saturation. 


In Eig. HI we show the normalized error, disturbance, 
LHS of the Branciard EDR in Eq. (flTl) . and LHS of the 
Heisenberg inequality in Eq. (TT^ . There are three char¬ 
acteristic regimes of the inequality, depending on the di¬ 
mensionless parameter a = Iquio /which deter¬ 
mines the intensity of a laser power or the strength of 
disturbance. At smaller cr, 5^ dominates the inequality, 
while at larger cr, both and Sr, increase in the same 
way. Since 5'e and Sr, are in the tradeoff relation for mid¬ 
dle O’, there exits the optimal a to minimise the LHS of 
the Branciard inequality at a certain frequency. One can 
show that the LHS is the closest to the saturation at the 
resonant frequency. The optimal disturbance strength 
there is given by CTopt = p-y/1 -b p^/16/(8-\/2) where 
p = KjuJm- With the optimal disturbance, the LHS of the 
inequality gives \/2 so that the inequality always holds 
but is not saturated. While the Heisenberg inequality is 
violated at around the resonant frequency in the small ct 
regime. This is understood from the inequality obtained 


by using the HUR, 5'e -b 5'^ > 2^^S^Sr, > 2 |x(D)|. For a 
damped harmonic oscillator, the lower bound is 2, while 
the minimum of LHS is S,: + Sr, = V2, which clearly 
violates the HUR. 

The gap ■\/2 from the saturation in the Branciard EDR 
arises from the fact that we have a factor of 2 in front 
of the Spg term in Eq. m- This is because of the dis¬ 
turbing force of the probe, Spg, included in both the er¬ 
ror and disturbance in Eqs. (l24l) and (l25l) . Since there 
are two tradeoff relations between Sgg and Spg (standard 
quantum limit 18, 1^) in the quantity S^ + S,,, both con¬ 
tribute to the lowest value of the LHS. One might think 
that the correlation between the probe variables Qq and 
i}), for example, squeezed light, helps reach the satura¬ 
tion. However, this does not work as we show above. It 
would be interesting to study whether this deviation from 
the saturation is specific to a damped harmonic oscillator 
or universal by applying to other systems. 

Conclusions — We have extended a universally valid 
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EDRs proposed for a discrete measurement to a contin¬ 
uous measurement and derived the EDRs in the Fourier 
space in terms of the power spectra of the system and 
probe variables. Application to a linear optomechani¬ 
cal system, particularly to a damped harmonic oscilla¬ 
tor, clearly shows the tradeoff relation between the error 
and the disturbance and leads to the existence of the 
optimal strength of the disturbance. The EDR in the 
Fourier space would be experimentally testable with, for 
instance, an optomechanical system in which zero-point 
fluctuations of an system contribute and would give a 
new insight into understanding of quantum foundations. 
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Supplemental material for ’’Universally Valid At the last line, we used 

Error-Disturbance Relations in Continuous 

Measurements” = A(—£7) , 


(38) 


Definitions of power spectra 


For an arbitrary observable A, the Fourier transform 
is defined by 

dO 

—. (33) 

-oo 

Then the one-sided power spectral densities are defined 
by 

(A(f7)Al'(f7') -b i:1'(f7')^(^^)) = 27r(5(f7 - n')Sx{n) , 

(34) 

(A(f7)y^(f7') -b Y'<{n')x{n)) = 27rS(fi - n')SxY(fi), 

which satisfy the properties 

Sxm = s*^(n) = Sx(-n), 

Srxm = S*xy{^) = SxY{-n) . 


If a system is stationary, the correlation function depends 
only on the difference of times t and t'. Denoting it by 
T = t — t', the correlation functions are related to the 
power spectral densities by the Wiener-Khintchine theo¬ 
rem: 

1 - - * ^ 1 dn 

-(A(r)A(O) + A(0)A(r)) = ^J ^ —e-*^^5y(D) , 

(35) 

1 1 do 

-(A(r)y(0) + r(0)A(r)) = - j ^—e-^^^SxY^ . 


In particular, taking r —^ 0 in Eq. (I35L we obtain the 
autocorrelation function 


{X\t)) 




dVt 

27r 


5jf(D) . 


Useful formulas 

We derive convenient formulas in the spectral domain. 
Formula 1: If X{t) is an Hermitian operator, then 

[A(D),At(D)] =0 . (36) 

proof 

X‘^{t)= [ — [ — . 

J 2tt J 2tt 

Integrating the above equation in time gives 

J dtx^{t) = J ^ J ^ x{n)x{n')2Trd{n + n') 

= J ^x{n)x{-n) 

= J ^x{n)xHn). (37) 


which is derived from Eq. (I33|) for an Hermitian operator. 
If we represent Eq. dSil) as the integral with respect to 
D', the expression becomes 

J ^vt( 0 ')A(D'). 

Since the variable X is an arbitrary at each frequency, 
we obtain 


A(D)A1'(D) = At(D)A(D) , 

that is, 

[A(D),a1'(D)] = 0 . 


Formula 2: For an arbitrary Hermitian operators, 
X{t) and Y{t), they satisfy the following relation in the 
Fourier domain, 

[Ao(H),yot(D')] = -2TrihSin-n'){RxYm - i?^x(D')} , 
with the response function defined in Eq. (15611 . 


proof 


[Ao(H),yo(f2')] 


dt 


' —OO J —OO 


dt' [Xo{t),Yo{t')] 


/ poo poo poo poo 

= f / dt dt'<d{t — t')+ / dt' / dtQ{t' — t) 

\J —OO J —oo J —oo J —oo 

xe“e*^'‘'[Ao(t),Fo(t')] 

= -ih{ f dt' RxY{n) 

L d — oo 

/ oo '' 

dte^^^+^''>*RYx{^') 

■OO y 

= — 27rihS(^D. + 12 ^) • 

Changing il' —)■ —il' and using Eqs. (t38|) and (113 , we 
have 


[Ao(D),yJ'(H')] = -2TTihd{n-n') {RxYi^) - R*Yxi^')} ■ 


Derivation of error-disturbance relations in 
continuous measurements 

We denote the system’s observables by A and B and 
suppose that the system’s variable A is measured by a 
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meter variable of a probe M. Following the Ozawa’s 
definitions in a discrete measurement @,01: we define 
two operators N and D in a continuous measurement: 


N{t) = Mit)-Ao{t), (39) 

D{t) = B{t) - Bo{t) , (40) 

where the subscript stands for free evolution. Then the 
error and disterbance induced by the measurements are 

el{t)^{N^t)) = {{Mit)-Ao{t)}^), (41) 

= (42) 


Here (• • •) denotes the ensemble average. The variances 
of observables Hq and Bq are 

alit) ^ {{Ao{t) - (io(t))}^) = {Alit)) - {Ao{t)r , 
al{t) ^ {{B,{t) - {BAt))?) = {Blit)) - iUt)? ■ 

Following the same proof of the universally valid error- 
disturbance relation by Ozawa in a discrete measurement 
@ 0 | , we show that the following error-disturbance rela¬ 
tion holds in a continuous measurement: 


iA{t)riB{t)+eA{t)<yB[t)+(TA[t)riB{t) > \CAB[t,t)\ , (43) 


where 

CAB{t.t')=^-{\AAt).m')]) ■ (44) 

proof 

Since M{t) and B{t) reside in different Hilbert spaces, 
then [M{t),B{t)\ = 0 gives 

[N{t), D{t)] + [N{t), Bo{t)] + [Mt),D{t)] 

= - [Aoit),Bo{t)] . 


Taking the expectation values and using the triangle in¬ 
equality, we have 


([iV(t),I)(t)]) + ([iV(t),4(t)]) + {[Ao{t),D{t)]) 


> 


([Ho(t), Bo{t)]) 


(45) 


From the definitions of the error and disturbance in 
Eqs. m and (HU), the expextation value of a squared 
operator is always larger than its variance. Then 


eA{t)>aN{t), riBit) > croit) . (46) 

Using the Robertson’s inequality 0, 


ax{t)(TY{t) > - ([Xo(t),Fo(i)]) 


which holds for an arbitrary pair of observables, we have 

eA{t)TjB{t) > i {[N{t),D{t)]) 


Substituting this for the first term in Eq. (I45|) and using 
the Robertson’s inequality and Eq. (gSl) again for the 
second and third terms in Eq. (14511 . we obtain 

eA{t)vB(t) + eA{t)cJB{t) + aA{t)r]B{t) 

>l\{[Ao{t),Bom)\- 


In a continuous measurement, it is convenient to work 
in the Fourier domain. To derive the error-disturbance 
relation in the Fourier domain, we assume that the whole 
system is stationary. For the later use, here we introduce 
a filter function T{t) and redefine the error and distur¬ 
bance in Eqs. (15^1) and (HOI) as 

Ar{t)= f dt'rit-t')N{t') = M{t)-Aoit), 

J —oo 

V{t)= f dt'r{t-t')D{t')=B{t)-BQ{t) , 

J —OO 


where Ao{t), Bo{t), M.{t), and B{t) are defined in a sim¬ 
ilar fashion using the same filter. The new error can be 
computed in the same way as the nonfiltered case. 


^ iJ^it)) 


/ dt' dt"T{t-t')r{t-t''){N{t')N{t")) 

J —OO J —oo 


where we defined 

/ OO pOO 

dr0(r)r(T)e*^^ = / drr(T)e*^^, 

-OO J 0 

and used Eqs. ((Ml) and dMO- As well, riB, 
are given by 

<^A = l ^|r(u)|^^.,(u), 

1 rIO 

4 = 2 j 


For the filtered quantities, the Ozawa’s inequality in 
Eq. (H51) is replaced with 

f-AlB + e-A(^B + (^AlB > \Cab\ ■ (47) 

Note that the time dependence in the original inequality 
disappears because we are assuming a stationary system. 















The right-hand side is 

\CAB\ = l\{[Aoit),Bom 


f dt' f dt"T{t - t')r{t -1") 

—oo —oo 




_ 1 

“ 2 

X e 




— oo 

—z(n^—o 




(48) 


If the system observables satisfy the commutation rela¬ 
tion with the following form 


[io(f^),-B^(o')] = 2TTih6{n - n')x{n), 


Eq. (H51) gives 


\Cab\ 


h 

2 



dVL' 

27r 


(x(o'))|r(fi')P 


(49) 


Substituting the above quantities for Eq. (|T7|) , we can 
obtain the Ozawa’s inequality in the Fourier space. How¬ 
ever, the expression includes many integrals and is dif¬ 
ficult to be dealt with. Since the band-pass filter we 
introduced can be chosen arbitrarily, we apply a narrow 
band-pass filter so as to extract the contribution from a 
certain frequency range centered at H with its 

width AH. Then the both sides of Eq. 621) gives the same 
factor |r(H)pAH/(27r), which are canceled out. There¬ 
fore, we finally obtain 


> ^l(x(f^))l • 

This procedure of the derivation can also be applied 
to other inequalities. We do not repeat the derivation, 
but we can express the Branciard’s inequality Q and 
the Robertson’s inequality [l^ in terms of power spectra 
as well. As a result, for a continuous measurement we 
obtain in the Fourier space the Branciard inequality 

+ 2^s,in)Sr,{n) {s.A^)s.Bm - l(x(f^))l'} 

>/^^l(x(^^))l^ ■ 


X{t) = Mt) + T / dt'[Hi{t')Joit)] + T 


and the Robertson’s inequality 


V'^.,(H)5.,(H)>ft|(x(H))| . 


A linear system 


If a system has Hamiltonian that is at most quadratic 
in its canonical coordinates and momenta, the system is 
called a linear system. Any linear combination of the 
canonical coordinates and momenta and a complex num¬ 
ber (c-number) defines a linear observable of the system. 
In the linear system, it is shown that in the Heisenberg 
picture, the commutator of the operators of any two lin¬ 
ear observables at two times gives a c-number [^. In 
what follows, we summalize the basic propaties of a lin¬ 
ear system. 

Suppose that A is an arbitrary observable of a system 
and that M is an arbitrary observable of a probe. We 
assume that the interaction between the probe and the 
system is the von Neumann type, whose Hamiltonian in 
the interaction picture is given by 


Hi{t) = -xo{t)Fo{t) , 


where xo{t) is the observable of the system to be mea¬ 
sured and Eo(0 is the generalized force of the probe. 
Since the variables in different Hilbert spaces always 
commute, the commutation relations between Aq (Mq), 
which is the freely evolving part of A (M), and the in¬ 
teraction Hamiltonian become 


[Hiit'),Ao{t)] = -2%C^A[t' , (50) 

[Hi{t'),Mo{t)\ = 2iCFM{t',t)xo{t') . (51) 


Here the correlation functions for freely evolving vari¬ 
ables are that defined in Eq. (l44l) . The time evolution 
of an arbitrary linear observable X{t) in the Heisenberg 
picture is given by, e.g. [^ . 





(52) 


In Eq. (I52|) . replacing X{t) with A{t) and M{t) and 


where Ao(t) is a freely evolving observable. 

















9 


using Eqs. dSp and flSlI) interatively, we obtain the time 
evolution of A and M 

A(t) = Ao^t) J dt'CAxit^t')F{t') , 

M{t) = Moit)+ % [ dt'CMF{t,t')x{t') . 

where the first terms are the variables in the free evo¬ 
lution and the second terms are contribution added by 
measurement procedures. 

When the above formulas are applied to a linear op¬ 
tomechanical system, in which the position of a system is 
measured by a optical probe, we identify the position and 
momentum of the test mass as A = and the position 
of the probe as M = Q. Then the dynamic of the sys¬ 
tem and the probe is governed by a set of the following 
equations: 

x{t) = Xo{t) + ^ / dt'Cxx{t,t')F{t') , (53) 

pit) =Poit) + dt'Cpxit,t')Fit') , (54) 

Qit) = Qoit) + T / dt'CQFit,t')xit') . (55) 

The Fourier transforms of Eqs. (Esii-dnsi) are 
i(f 2 ) = io(f^) + Rxxi^)FiVL) , 
m) = poi^) + RpxAt)F{n) , 

Q(o) = QoAt) + i?Qi7’(n)i(n), 


with quadrature annihilation operators 

_ a+ -I- oL _ a+ — oL 

= -TT- • ■ 

where A is an effective beam area, wq is the carrier fre¬ 
quency, and c is the speed of light. The annihilation 
operators, aj- = a(a;o ± fl), satisfy the commutation re¬ 
lations 

[a±{n), aiin')] = 2TTSin - n '), 
[a±(f7),a±(0')] = 0 , 

then the quadrature annihilation operators satisfy 
[oi(17), 4(17')] =-[a2(17),4(f^')] = 2f7r(5(17-17') , 
[ai(17),ai(17')] = [m(17), 4(17')] = [ai(17), 02(17')] = 0 . 

We suppose that light is reflected at a test mass, e.g. 
a tiny mirror whose size is neglected. We denote the 
small displacement of the mirror by x = cAt/2 (here 
we treat the mirror dispalcement as a c-number only for 
the purpose to derive the relevant quantities of the light 
probe). Writing an input field as a and a output field 
as b and substituting Eqs. (|5^ and (|5^ into the rela¬ 
tion E},{t) = Ea{t — At) with the assumption that the 
sideband frequency is much smaller than that of a car¬ 
rier, that is, cjQ ^ 17, we have the input-output relation 
for the quadrature operators at the leading order in the 
small quantities 17/w, x and 01 . 2 , 

61 = Oi , 

62 = 02 -f , (60) 


where we defined 


2 /■“ 

7?xv(17) = - / dr 0(r)C'xv(r)f 

^ J —00 

poo 

dr CxY{x)e'^^^ . 


,iVtT 


_ 2 
h 

From the definition of Rxvi^), we have 

Rxvi^) = 7?xy(—17) . 

Optical probe 


(56) 


(57) 


Suppose that carrier light with the classical amplitude 
Eq and the frequency ujq is injected to measure the posi¬ 
tion of an object. According to the two-photon formal¬ 
ism in 2 l|, an electromagnetic field of the carrier light 


including vacuum fluctuations is written as 


Eait) = [Eo+Eaiit)] coswof + Ao 2 (t) sinwot (58) 

dQ 

’ 

(59) 


Eal.lit) — 


Ac 


(ai ,2 e 


— iQ,t 


2 6 


iClt 


where the DC component of the injected laser power is 

EnAc 


do = 


Stt 


From Eq. (l60l) . we identify the position operator for the 
electromagnetic field as 


Qoidt) — 


' hc^ 
SIqujo 


02 ( 17 ) . 


(61) 


Next let us calculate radiation pressure exerted on the 
test mass. The radiation pressure is given by 

4p(t) = ^ = Asiit) , 

C ZTT 

where I is an optical power injected to the test mass. 
The factor 2 comes from that the light is reflected by 
the test mass. By smoothing out the rapidly oscillating 
terms with the frequency wq, the fluctuating component 
with the period of a sideband is 

A 


AF^pit) — —EoEaiit) 
Ztt 


SIofkjjQ 


oie 




27r 
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Thus we define the stochastic force exerted by the mea¬ 
surement 

Fo{n) = . (62) 

From Eqs. (El]) and (162]) . we have 

[Qo{n),F^ (11')] = -^TTihSin - n '), (63) 

[go(f^),QS(f^')] = [Fom,F^m = o- 

Since the Fourier transform of Eq. (l63ll is 

[go(t),^o(i')] =-ihS{t-t') , 

by writing r = t — t', the response function in the Eourier 
space is 

RQFin) =-j dTQ{T)CQF{T)F^^ 

= 1 . (64) 

Using Eq. (IMl) and assuming the sidebands of the 
probe is in a coherent vacuum state (no correlation be¬ 
tween Q and F), namely, 

(ai(fl)a](ll') -I- a](H')ai(U)) = 2f5{Q, — , 

we obtain the power spectra for the coherent light probe 

= SFom = ^^, SF,Q,m = o. 

OlQiOo c 


Braginsky’s inequality 

In a linear system, all commutation relation gives a c- 
number. Suppose that the Hermitian variables A and B 
of the probe satisfy the commutation relations: 

[i(u),.Bi'(u')] = 2 TTihx{fi)Sin - n'). ( 65 ) 

The Braginsky’s inequality can be derived from the com¬ 
mutation relation in Eqs. (1651) and the Schwarz’s inequal¬ 
ity for anbitrary quantum states, |a) and \f3): 

{a\a){m > \{am^ . 

Let us consider the case with ja) = A(U)|-) and 
1/3) = .B(fl)|-), where I-) denotes an arbitrary state. The 
Schwarz’s inequality gives 

{■\AHn)Am){-\B^in)Bm-) > |(.|it(ii)F(u)|.)p. 

( 66 ) 

From Eq. (1361) . the Hermite-conjugate operators in the 
Fourier space commute, 

(.|it(U)i(U)|.) = i(.|i(ll)it(U) +it(ii)i(U)|.) 

= i2^5(0)FA(L!) , 

(•|#(U)F(U)|.) = i(.|B(L!).Bt( 0 ) + #(U)F( 0 )|.) 

= ^2^5(0)55(11) , 


On the other hand, 

{■\AAn)Mm-) 

= i(.|it(ll).B(ll) + .B(ll)it(ll) + [it(o), B(ll)]|.) 

= l 2 nS{ 0 )SBA{n) + l{-\[AA^),Mm\-) 

= i2^5(0)55A(f^)-i(-|[i(f^),Bt(ll)]t|.) 

= ^27rS{0) {SBAm + tnx*m} , 

where we used Eqs. (l65|) . Substituting these into 
Eq. (1551) . we obtain 

5a(11)5b(11) 

> |55A(11) + ^%*( 11 )|" 

= |55a(^^)P + + 2?iIm[x(ll)55A(^^)] . 

Thus, we have 

SAmsBm-\SBAm^ > h^\xm^+‘^fiMx{^)SBAm. 

(67) 


Derivation of Eq. (1321) 


The LHS of the Branciard inequality in the case of a 
damped harmonic oscillator is given by 

Q I c _ ^Qo F ‘^\Rxx\’^S fq + 2 Re[i7a;3;55QQo] 

' 2ftIm[J?,J - • 

We assume that SfqQo is real-valued. Using the Bragin¬ 
sky’s inequality Sq^^Sfo — in Eq. (l67l) . we 

have 


Se + Sr, > 


2^2\R,,\^Sq„Sfo + 2Re[Jt;,,^5F„Qj 
2ftlm[i?3,x] 


2v'2|i7,;,3;p(fi^ + |5 FoQoP) + 2Re[Fa;3;5FQQp] 
2hlm[Rxx] 


= ^/2(1 -I- a2)(l -I- 62) -I- ab 
> \/2(l -I- a2)(l -I- 5^) — Va?b‘^ 


where 


Since 


^_ SfQq ^_ lLe[l?a, 3 ,] 

h ' ~ lm[Rxx] ' 


(\/2(l + a2)(l + 62) _ v^)2 

= 2(1 -f aAA + bA + a^b^ “ 2y/2a^b^{l + a^){l + b^ 
> 2{1 + aAA + bA + a'^b^ 

-V2{a^{l + bA + bAl + aA} 

= 2 +{2- V2){a^ + bA + (3 - 2V2)a^b^ , 



















11 


the LHS of the Branciard inequality is constrained from 
below 


S, + Sr, > y^ 2 + (2 - \/2)(a2 + 62) + (3 - 2V2)a%^ 

> V2 . 



